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An  expression  for  the  discontinuity  of  theexchange  cor¬ 
relation  potential  v  of  an  Insulator  is  derived  entirely 
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within  the  framework  of  density  functional  theory.'^The  dis¬ 
continuity  is  expressed  in  terms  of  changes  of  the  exchange 
correlation  energy,  E  ; ,  of  a  perfect  N-particle  insulator 
when  (a)  a  conduction  electron  is  introduced  (b)  a  valence 
electron  is  removed  and  (c)  -Wm  external  perturbation  is 
applied  to  the  perfect  insulator  (without  changing  N)  such 
tnat  the  density  change  is  equal  to  minus  the  sum  of  the 
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Since  the  work  of  Perdew  and  Levy'  and  of  Shan  and 
Schliiter*  it  has  been  known  that  the  exchange  correlation  po¬ 


tential  Kohn-Sham  (KS)  equations’  has  an  r-in- 

dependent  discontinuity  A,  as  one  crosses  the  energy  gap. This 
discontinuity  has  been  studied  with  the  aid  of  Green's  func¬ 
tion  theory  by  Wang  and  Pickett**,  Sham  and  coworkers*  »  *  '' 
and  Hanke*.  In  this  note  we  discuss  the  discontinuity  enti¬ 
rely  from  a  density  functional  view  point. 

We  shall  consider  an  Insulator  at  temperature  T«0‘*’K  in 
three  physical  ground-states:  1.  with  N  electrons,  without 
conduction  electrons  or  holes.  We  denote  the  corresponding 
density  distributions  by  no(r)  with 

/n  0  (r  )  dr  -  N  ;  (  1  ) 


2.  with  N+v  electrons;  and  3.  with  N-v  electrons,  with  v<<N.’ 
We  write  the  densities  in  states  2.  and  3.  as 

n  ■*■  (  r  )=n  0  (  r  )  +  V  n'^(r);  /n'^(r)dr  »  1. 

(2) 

V  V 

n7'!"'=no;r)-v  n  (r);  /n  (r)dr  »  1. 


We  take  (kT)  much  smaller  than  any  physical  energy  but  larger 
than  the  energy  spacing  between  successive  single  particle 
excitations.^  The  system  is  in  contact  with  a  particle  bath, 
allrwlng  continuous  changes  of  the  total  particle  number'. 

5y  its  definition,  the  energy  gap  can  be  expressed  in 
t-^rms  t  ground  state  energies. 


Each  of  these  energy  is  given.  In  the  KS  theory’,  by  the 
express  ion 


Ev  [n  (r  )  ] -T  [  n  ( r  )  ]  + /v  (r  )n  (r  )dr  *  ^  f  ---7-777" 

V  S  j  *  *  I 


♦E^^[n(r )  ] , 


(it) 


where  v(r)  is  the  (fixed)  external  potential  and  the  other 
symbols  have  their  asiial  meanings.  As  is  well-known,  this 
expression  can  be  transformed  into 


Ey  [n(r  )  ]-[e  Jn(r)  J*  ^^7x77-- dr  dr  ’  +  [n  ( r  )  ] 


(5) 


-/  V ( r ' ; [ n (r )  ]  )n ( r ’ ) dr ’ 


where  i  runs  over  al-1  occupied  single  particle  levels  associ 
ated  with  the  KS  equation 

t-1^  V^+v(r)+  /  “£7"  dr  ' +v  (r  ;  n  [r  ’  ] )  }  ^  (  r  )  -  e  ^  ^  (  r  )  ,  (6) 

and 

V  (r;[n(r')])=6E  [n(r’)]/5n(r).  (7) 

xc  xc 


We  have  emphasized 
‘ i 0 h a  1 s  of  n ( r  '  i  . 


in  our  notation  that  v  and  e, 

xc  i 


are  func- 


The  difference  E,,  -E.,  in  (3)  can  now  be  calcilated  from 

N  +  v  N 

(5),  to  first  order  In  v: 


N  +  v  N 

^N  +  v7^N  ■  ^  ^i  J  ’^‘^i  “ 

N  +  1  I 


n  C r  )n  0  ( r  '  ) 


dr  dr  ’ 


♦  6Ej^^  [n  (r  )  ]“  V  /  v ( r '  ;  [ n o  ( r  )  ]n ( r  '  ) dr  '  - 


/  6v ( r ' ; [n (r )  ]  )n 0 ( r  ’  )  dr’ 


Here  Is  the  lowest  conduction  band  energy  of  the  KS 

equation  for  the  N+v  particle  system  (v+O),  calculated  with 
the  exchange  correlation  potential  ^^^(r)  appropriate  for 
this  system; 

c  6E  [no(r)+vn®(r)] 

v""  (r)  =  lim  — - t-T-T -  (9) 

xc  6n(r) 

v+0 


The  simplification  in  the  last  step  of  Eq.  (8)  is  due  to  the 

cancellation  of  the  2*^^,  and  last  terms,  and  of  the 

t  h 

and  5  terns.  Similarly 


'n'^N-v 


'N-l 


where  is  computed  with 


■5  E..  _  [  n  0  (  c  )*■  vn'^  (  r  )  ] 


Thus  the  gap  Is  given  by 


E 

g 


c 


e 


V 

N-1  * 


''12) 


The  arguments  of  refs.  (1)  and  (2)  show,  and  we  shall  verify 
that 


a  constant  independent  of  r.  Therefore  (12)  can  also  be  writ' 
ten  as 


•N>1 


’N  +  l 


-  e 


N**1 


^g 


(  U) 


where  e  is  the  non*physical  gap  of  the  KS  single  particle 
g 

V  c 

Insulator,  computed  with  either  v  or  v  , 

xc  xc 

Q 

We  now  turn  to  a  consideration  of  Eq.  (9).  We 

introduce 


6n^,(r)  =  y  5(r*r')  ,  /  6(r-r')  dr  ••  1  (15) 


where  Y  <<  1  and  5  is  a  normalized,  regularized  6*"function 


(r')  -  lim  lira  ^ I [ n ^ ( r )  +  v  n  ( r  )  + Y6 ( r-r ' ) ]  - 


v-^o  y*Q 


“^xc  ] } 

is  now  useful  to  decompose  6(r“r*)  into  two  parts:  n 
which  increases  the  number  of  electrons  by  1 ,  and  a 

remainder,  m^,(r),  which  leaves  the  number  of  electrons 
unchanged  at  N. 


5(r-r')  «  n  (r)  +  m^,(r), 


where,  evidently  in  view  of  Eqs.  (2)  and  (15), 


/  m^,(r)dr  =  /[6(r-r')  -  n°(r)]d 


r  -  0. 


Substituting  (1?)  into  (16)  gives  two  terms. 


V  (r  '  ) 
xc 


c  c 

y  +  w^  ( r '  )  , 

XC  X  c 


where 


^xc  =  -  E  [no(r)] 


'  X  c  ,  M  +  1  X  c  ,  N 
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and 


w 


0 

xc 


(r*  ) 


lim  lim  |  E  [  n  « ( r  )  +  vn°  (  r  )  +  Ym°  ( r  )  ]-E  [  n  o  (  r  ) 

\;->0  Y->0  '  r  xc 


+  vn 


(r)]l 


(21  ) 


The  two  density  arguments  in  Eq .  (21)  differ  by  Y  m^,(r), 
defined  by  Eq.  (17),  and  corresponding  to  6N«0.  This  density 
difference  must  therefore  be  understood  as  being  brought 
about  by  the  action  of  a  small  external  perturbing  potential, 
Y  u^,(r)  modifying  the  (N+v)  particle  ground  state.  In  the 
limit  v->0  the  role  of  the  conduction  electrons  in  (21)  beco¬ 
mes  negligible,  so  that 


xc 


(r')  -  lim  •7(E^^[no(r)  +  Ym^,  (r  )  ]-E^„[no(r)  ]  }  .  (22) 


Y-»0 


xc ' 


In  a  completely  analogous  manner  we  obtain  the  following 

V 

resul ts  for  v  ( r ' ) : 

xc 


(  r  '  )  = 


xc 


W  (  r  ’  ) 
xc 


(23) 


where 


xc 


E,,ln,(r!l  -  E^Jn.(r)  -  n^r)]  , 


'xc  ,  N  "  ^xc  ,N-  1  ’ 


(211) 


-  8  - 


r 


(r”)  -  lim  YU„-[no(r)]  -  E  [no(r)  -  ,  ( r  )  ]  (25) 

xc  Y“^0  *  ^ 


with 


m^,(r)  =  6(r-r*)  -  n'^(r). 


From  the  expressions  (19),  (20),  (22)  and  (23),  (2U), 

C  V 

(25)  we  can  calculate  the  difference,  v  -v  .  Note  that  the 

XC  ‘  xc 

two  particle  number  cons erv ing  changes,  (22)  and  (25),  can  be 


combined 


^xc^'^“^’'^*''^'"r'  ^'"^:"'r’  ]"Ej^^,[no(r)  ] 


independent  of  r'.  Thus  we  obtain  for  the  discontinuity  of 


't  > 


“  9  - 


’.‘•Vv' 


Exc["o(i'  )+n°  (r)  ]-Ex^[n«(r  )  ]  }  +  lE^p[no  (r  )-n''  (r  )  ]-E^^  ]  } 


+  |Exp[no(r)-(n°(r)-n''(r))]-Exj,[no(r)}.  (28) 


These  three  terms  represents  the  changes  in  the  exchan¬ 
ge-correlation  energy  of  the  N-partlcle  Insulator  when  (a)  a 
conduction  electron  is  added;  (b)  a  valence  electron  is  remo¬ 
ved;  and  (c)  an  external  potential  is  applied  which  changes 
the  density  of  the  N-particle  insulator  by  ■*•  ( n  (r)^n  (r)) 
without  introducing  either  electrons  or  holes.  Note  that  if 
Exc[n(r)]  had  a  regular  dependence  on  n(r)  near  n(r)«no(r), 
then  the  three  terms  could  be  expanded  in  the  small  quant i- 

C  V 

ties  n  (r)  and  n  (r)  resulting  in 


oE  [n(r')'l 
xc  '■ 


r eg  on ( r  ) 


n ( r  )  »n 


{n^(r)-n'^(r)-(n^(r)n^(r))  }dr-0 


Thus  any  approximate  theory  which  uses  a  regular  expression 
for  ^xc”  ^ example,  the  local  density  approximation  - 

must  yield  a  vanishing  A.  The  correct  formal  expression  (28) 
makes  it  clear  however,  on  physical  grounds,  why  A  does  not 
vanish:  the  effect  on  E^^  of  modifying  the  external  potential 
of  the  insulator  so  that  its  density  changes  by 


v.v. 


V/.vy, 


-  10  - 


C  V 

-(n  (r)-n  (r))  does  not  cancel  the  sum  of  the  physically 
totally  different  changes  of  adding  an  electron  with  density 

C  V 

n  (r)  and  of  removing  an  electron  with  density  n  (r).  The 

correct  E^^[n(r)]  will  give  the  correct  A  by  Eq.(28).  It 

remains  a  challenge  to  find  useful  non-regular  expressions 

for  E  which  will  yield  accurate  values  for  A. 
xc 
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